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The forces on the vortex, transverse to its velocity, are considered. In addition to the superfluid 
Magnus force from the condensate (superfluid component) , there are transverse forces from thermal 
quasiparticles and external fields violating the translational invariance. The forces between quasi- 
particles and the vortex originate from interference of quasiparticles with trajectories on the left and 
on the right from the vortex like similar forces for electrons interacting with the thin magnetic-flux 
\Q • tube (the Aharonov-Bohm effect). These forces are derived in the Born approximation for phonons 

OA ' from the equations of superfluid hydrodynamics, and for BCS quasiparticles from the Bogolyubov-de 

ON ■ Gennes equations. 

The effect of external fields violating translational invariance is analyzed for vortices in the two- 
dimensional Josephson junction array. The symmetry analysis of the classical equations for the 
array shows that the total transverse force on the vortex vanishes. Therefore the Hall effect which is 
linear in the transverse force is absent also. This means that the Magnus force from the superfluid 
component exactly cancels with the transverse force from the external fields. 
The results of other approaches are also brought together for discussion. 

> 1 i. introduction 

On 

The Magnus force on a vortex has long been known in classical hydrodynamics . This force appears if the vortex 
\Q . moves with respect to a liquid. The force is normal to the relative vortex velocity and therefore is reactive and does 
not produce a work. In general, such a force arises always when a body with a flow circulation around it moves 
through a liquid or a gas (the Kutta- Joukowski theorem) . The most important example is the lifting force on a wing 
of a aeroplane which keeps the aeroplane in the air. 

The key role of the Magnus force in vortex dynamics has become clear from the very beginning of studying superfluid 
hydrodynamics |^||. The superfluid Magnus force was defined as a force between a vortex and a superfluid and 
therefore was proportional to the superfluid density p s . But in the two-fluid hydrodynamics the superfluid Magnus 
force is not the only force on the vortex transverse to its velocity: there was also a transverse force between the 
vortex and quasiparticles moving with respect to the vortex. The transverse force from rotons was found by Lifshitz 
and Pitaevskii j| from the quasiclassical scattering theory. Later Iordanskii || revealed the transverse force from 
' phonons which was equal in magnitude and opposite in sign with the quasiclassical force of Lifshitz and Pitaevskii. 
, From the very beginning the Iordanskii force was a controversial matter. Iordanskii suggested that his force and the 
Lifshitz-Pitaevskii force were of different origins and for rotons they should be summed. As a result, he concluded that 
the transverse force from rotons vanished. But the analysis done in Ref. j(| demonstrated that the Iordanskii force 
for rotons is identical to the Lifshitz-Pitaevskii force and they must not be added. In addition, the Lifshitz-Pitaevskii 
force from rotons was calculated in the original paper Q with a wrong sign. After its correction the transverse 
force on the vortex had a same sign and a value both for rotons (the Lifshitz-Pitaevskii force) and for phonons (the 
Iordanskii force). In the same paper Q it was pointed out that the Iordanskii force for phonons and rotons results 
from interference between quasiparticles which move past the vortex on the left and on the right sides with different 
phase shifts, like in the Aharanov-Bohm effect 0. 

In the theory of superconductivity the Magnus force appeared first in the paper by Nozieres and Vinen j^]. In 
clean superconductors the BCS quasiparticles produce an additional transverse force on the vortex P,|lO| analogous 
to the Iordanskii force in superfiuids. The total transverse force is responsible for the Hall effect in the mixed 
state of a superconductor. But the Hall effect was rather weak in classical superconductors. An explanation of it 
was suggested by Kopnin and Kravtsov [|llj (see also Ref. [Q): impurities interact with quasiparticles bound in 
the vortex core and this interaction produces an additional transverse force on the vortex. In contrast with the 
quasiparticle transverse force which increases the total transverse force, the impurity force decreases it and in a dirty 
superconductor completely cancels the Magnus force. As a result, the strong Hall effect is possible only in superclean 
superconductors. The transverse force from the bound states in the core has been recently rephrased in terms of the 
spectral flow through the quasiparticle bound states fl3}| . 
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A new wave of interest to the Magnus force came with discovery of high-T c superconductivity. A few reasons of this 
interest might be mentioned: (i) The so-called Hall anomaly was observed JlJ]: near T c the sign of the Hall voltage 
is opposite to that expected from the standard vortex dynamics, (ii) It has become possible to obtain superclean 
single crystals with large Hall angle as predicted in the superclean limit of the theory of Kopnin and Kravtsov. This 
made possible to observe magnetoresonances in the a.c. response of the superconducting single crystals connected 
with waves propagating along vortices jl5|,[l6| . (iii) The effective Magnus force governs quantum vortex nucleation in 
clean high-T c superconductors intensively discussed now [l7|Jl^ |. 

Despite a lot of work done to understand and calculate the Magnus force, it still remains to be a matter of controversy 
with a number conflicting points of view on it. A new discussion has been launched by the paper of Ao and Thouless 
|jl9f . The main claim of Ao and Thouless is that there is an universal exact expression for the total transverse force on 
the vortex (the effective Magnus force) which does not depend on the presence of quasiparticlcs or impurities. This 
force derived from the concept of the geometrical phase (the Berry phase) coincides with the superfluid Magnus force 
and therefore is proportional to the superfluid density. According to Ao and Thouless, there is no transverse force 
on the vortex from quasiparticles and impurities, though they might influence the value of the superfluid density and 
thereby influence the amplitude of the Magnus force. 

The Ao-Thouless theory is in an evident disagreement with the previous calculations of the transverse force on the 
vortex (the effective Magnus force) in superfluids and superconductors reviewed above. It attracted agreat attention 
and has been supported in a number of recent publications of other theorists (see e.g. Refs. p0|-p2[). If the Ao- 
Thouless theory were true, it would be necessary to revise the whole basis of the vortex dynamics. For example, 
on the basis of the Ao-Thouless theory Simanek [EH] suggested that quantum vortex tunnelling is governed by the 
Magnus force obtained from the Berry phase approach, i.e., proportional to the superfluid density, in contradiction 
to the previous theory |l7],[ll|. Therefore it is important to understand what the Magnus force is and whether the 
Ao-Thouless theory is true or not. 

The present paper is to analyze this controversy. Among the sources of controversy there is semantics. Therefore 
it is important to define force terminology from the very beginning. The word force itself is only a label to describe a 
transfer of momentum between two objects. Before using these labels one must to analyze the momentum balance for 
any object and only then to label various contributions to these balances as forces. Keeping this in mind, the forces 
under discussions may be defined as the following: 

• There is a momentum transfer between a superfluid and a vortex. It is revealed analyzing the momentum balance 
for the superfluid moving with the velocity v s whereas the vortex moves through the superfluid with a different 
velocity vl- This momentum transfer is the superfluid Magnus force. It is proportional to the superfluid density 
p s and transverse to the relative velocity v*l — v s . 

• Analyzing the momentum balance for the whole translationally invariant liquid (including the superfluid and 
normal parts of it) around the vortex one may reveal a contribution presenting the momentum transfer between 
the vortex moving with the velocity vl and the normal fluid (the gas of quasiparticles) moving with the velocity 
v n . This force is proportional to the relative velocity vl — v n and has the components longitudinal and transverse 
to Vl — v n - The transverse component is the Iordanskii force . 

• If there is no translational invariance, as in a dirty superconductor, the momentum balance for the whole 
liquid must include also forces external for the liquid, namely, the momentum transfer to the impurities rigidly 
connected with the crystal lattice. When the latter is at rest, this momentum transfer, or the force from 
impurities, is proportional to the absolute value of the vortex velocity vl- Its component transverse to vl is the 
Kopnin-Kravtsov force. 

• The momentum balance for the whole liquid around the vortex is at the same time an equation from which 
one must find the vortex velocity Vl- Therefore it is useful to collect all the terms proportional to vl together. 
After it the term uniting all contributions transverse to vl is the total transverse force on the vortex, or the 
effective Magnus force. 

Thus in general three forces contribute to the effective Magnus force: the superfluid Magnus force, the Iordanskii 
force, and the Kopnin-Kravtsov force. Experimentally they can determine the effective Magnus force, but not the 
"bare" superfluid Magnus force. In rotating superfluids the effective Magnus force determines the mutual friction. 
One can find the latest reviews of the experiment and the theory on mutual friction in 3 He in Refs. [^3[^4| . In 
superconductors the Hall conductivity P~lJ2q] and the acoustic Faraday effect for the transverse ultrasound wave 



propagating along vortices 26 are linear in the effective Magnus force. The process of vortex quantum tunnelling is 
also influenced by the effective, but not the superfluid Magnus force. So the final outcome of the theory must be the 
amplitude of the effective Magnus force. Its presentation as a combination of three forces is an intermediate stage of 
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the theory. In fact, this presentation is valid only if (i) the number of quasiparticles is not too large and their mutual 
interaction is not strong; (ii) external fields violating translational invariance are not too strong. The first condition 
is violated close to T c , where other approaches based on the Ginzburg-Landau theory (or its analogue for superfluids, 
the Ginzburg-Pitaevskii theory) must be used [ ^7|j28| . The second condition doesn't hold in the Josephson junction 
array considered in the present paper (see below). In these cases the theory deals directly with the effective Magnus 
force in the equation of vortex motion: its decomposition on the "bare" superfiuid Magnus force and the forces from 
quasiparticles or impurities becomes conventional and of a little physical sense. 

Whereas there is a consensus among theorists on the superfiuid Magnus force, the Ao-Thoules theory rejects the 
Iordanskii force from quasiparticles and the Kopnin-Kravtsov force from impurities claiming that amplitudes of the 
effective and the superfiuid Magnus forces are exactly equal. Therefore the present paper considers the effect on 
the force balance from (i) quasiparticles, and (ii) from external fields violating the translational invariance of the 
superfiuid. 

In my analysis of the quasiparticle effect I chose the phonon-vortex interaction which may be described by the 
nonlinear Schrodinger equation long ago suggested for a weakly nonideal Bose-gas (the Gross-Pitaevskii theory Q ) ■ 
It is well known that the nonlinear Schrodinger equation yields the usual superfiuid hydrodynamics. It is a good 
starting point for further discussion, which one may expect a consensus of all parties on. Indeed, according to 
Aitchison, Ao, Thouless, and Zhu 30 1, the nonlinear Schrodinger equation, and the superfiuid hydrodynamics derived 
from it, is a reliable model to describe both Bose- and Fermi-superfluids near T = 0. The next step is to analyze 
scattering of the sound wave (phonon) by the vortex in hydrodynamics. Just at this stage a disagreement appears. Ao 
and Thouless believe that this scattering can produce only a dissipative force on the vortex, but not a transverse one. 
Recently Demircan, Ao, and Niu (3^] tried to prove it using the Born approximation. But they ignored peculiarities 
of the phonon Born scattering at small angles which resulted in the Iordanskii force. It is important to note that the 
controversy arises not from a difference in ideology: anyone is free to choose a language to derive the Magnus force 
at T = 0: either the standard hydrodynamics, or newest topological concepts of the geometrical phase. But there 
is a disagreement in calculation of integrals describing the phonon scattering in the first order of the perturbation 
theory. We hope to show in this paper at which point Demircan, Ao, and Niu [j32| missed to take into account the 
Aharonov-Bohm interference of phonons which was ignored by the Ao- Thouless theory. 

Now discussions around the transverse force on the vortex in the presence of impurities violating the translational 
invariance concentrate mostly on the contribution of the core bound states in Fermi superfluids. This requires a 
rather sophisticated analysis (see [fl3f and references therein). In the present paper I chose another example when 
the translational invariance is absent: the two-dimensional Josephson junction array(JJA). This is a regular lattice 
of nodes with the Josephson coupling between them. Experimentally, any node corresponds to a superconducting 
island in an artificially prepared JJA, or to a grain in a granular superconducting film. The behavior of the JJA 
in an external magnetic field is usually described in the picture of moving vortices similar to the mixed state of 
type II superconductors. The dynamics of vortices in JJA attracts a great interest of experimentalists [[33 34 and 



theorists |35| |4J]. There is an intrinsic pinning of vortices at the JJA cells, and vortices can move only if the driving 
supercurrent is more than the critical value. But when they start to move, in many cases a good approximation is 
to replace the lattice by a continuous superconducting film. However, the hydrodynamic derivation of the Magnus 
force is not valid since it assumes the momentum conservation law and the translational invariance. In the present 
paper it will be shown that the Hall effect is exactly absent in the classical theory of JJA which neglects the charge 
quantization. Since the Hall effect is linear in the amplitude of the effective Magnus force, the latter also vanishes in 
the classical JJA. This statement directly follows from the symmetry of the dynamic equations. At the same time the 
superfiuid density is finite in the continuum limit of JJA and therefore the superfiuid Magnus force doesn't vanish. 
Therefore the theory based on the Berry phase approach predicted a finite effective Magnus force and the Hall 
effect for JJA. Our result might be interpreted as that the superfiuid Magnus force is compensated by some force 
external for the liquid, like the Kopnin-Kravtsov force in a dirty superconductor. But as mentioned above, JJA is 
a system with a strong violation of translational invariance, for which this interpretation is purely formal. Only the 
resultant effective Magnus force has a physical meaning. 

We start from Sec. || which shows how the Magnus force appears in the phenomenological theory of neutral and 
charged superfluids. The force terminology is also introduced explaining to which term and in which equation any 
force under discussion corresponds. Section III is devoted to the transverse force bet ween quasiparticles and the vortex 
(the Iordanskii force) and its connection with the Aharonov-Bohm effect. Section III A| reminds connection between 
the nonlinear Schrodinger equation for the condensate (the Gross-Pitaevskii theory) and superfiuid hydro dynamics 
and phonons. Scattering of the sound wave (phonon) by the vortex in hydrodynamics is analyzed in Sec III B using the 
Born approximation and the effective cross-sections. It is shown that the standard scattering-theory approach fails to 
reveal the transverse Iordanskii force because of the divergence of the scattering amplitude at small angles of scattering. 



The analysis of the small-angle scattering is presented in Sec. Ill C . It reveals the interference between quasiparticles 



with trajectories on the left and on the right from the vortex. In Sec. Ill D the same results is rederived using the 
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partial-wave expansion, and the analogy with the Aharonov-Bohm effect is shown. A more general quasiclassical 



derivation of the transverse force from the quasiparticles with an arbitrary spectrum is presented in Sec. Ill E. In Sec 



IV the transverse force between the BCS quasiparticles and the vortex is derived using the Born approximation for 
the Bogolyubov-de Gennes equations. Section |Vj presents the symmetry analysis of the classical dynamical equations 
for JJA which shows that the total transverse force on the vortex in JJA vanishes and as a result of it the Hall effect 



is possible only in the quantum theory of JJA which takes into account charge quantization. The last section VI 
contains the summary and the discussion of other approaches to the problem. 

II. WHERE AND HOW THE MAGNUS FORCE APPEARS 

A. The Magnus force in classical hydrodynamics 

As mentioned in Introduction, the Ao-Thouless theory connects the Magnus force with the concept of the geometrical 
phase. But for better understanding of the Magnus force origin it is worth to remind how the Magnus force arises in 
classical hydrodynamics. 

Let us consider an isolated straight vortex line in an incompressible inviscid liquid. The line along the axis z induces 
the velocity field 

, _» ac x r 

Mr) = ^ • (1) 

Here r is the position vector in the plane xy, and k is the circulation vector directed along the axis z. The circulation, 
given by 

k = <j> v v ■ dl , (2) 

may have arbitrary values in classical hydrodynamics. In addition, there is a fluid current past the vortex line with a 
velocity vq. Then the net velocity field around the line is 

v(r) = v v (r) + v . (3) 

The Euler equation for the liquid is 

dv - 1 - F 

— + (v-V)v = —VP + -6 2 {7 f ). (4) 
at p p 

Here p is the liquid density and P is the pressure. This equation suggests that an external (5-function force F is 
applied to the liquid at the vortex line. 

Assuming that the vortex line moves with the constant velocity vl, one obtains that 

^ = -(vl-^)v. (5) 

Then the Euler equation (JI]) yields the Bernoulli law for the pressure: 

P = Pa - ^p[v{r) - v L } 2 =Pq - -pv v (r) 2 - pv v (r) ■ (v - v L ) . (6) 

Here Pq and Pq = Pq — ^p(vo — ^l) 2 are constants which are of no importance for the following derivation. 

Next one should consider the momentum balance for a cylindrical region of a radius tq around the vortex line. The 
momentum-flux tensor is given by [EBJ 

Hy = PSij + pVi(r)Vj(f) , (7) 
or in the reference frame moving with the vortex velocity vl- 

%j = P*ij + ~ v Li ){v 3 - v Lj ) . (8) 
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The momentum conservation law requires that the external force F on the vortex line must be equal to the 
momentum flux through the entire cylindrical boundary in the reference frame moving with the vortex velocity vl- 
The latter is given by the integral J dSjU'^ where dSj are the components of the vector dS directed along the outer 
normal to the boundary of the cylindrical region and equal to the elementary area of the boundary in magnitude. 
Then using Eqs. (0), @, and (||), the momentum balance yields the following relation: 

p[(v L - v ) x Ft] = F . (9) 

On the left-hand side of this equation one can see the Magnus force as it comes in the classical hydrodynamics. A 
half of this force is due to the Bernoulli contribution to the pressure, Eq. (jfj); another half is due to the convection 
term oc ViVj in the momentum flux. The Magnus force balances the resultant of all external forces applied to the 
liquid at the vortex line (the force F). In the absence of external forces the vortex moves with the velocity of the 
liquid: vl = Vq (the Helmholtz theorem). 

This derivation demonstrates the classical origin of the Magnus force: quantization of circulation is not necessary 
for its existence. During the derivation we referred to the hydrodynamic equations only at large distance ro S> r c 
from the vortex line. It might seem as if the fluid in the vortex core did not matter at all. However, the derivation 
is based on the assumption that the momentum is a well-defined conserved quantity everywhere even inside of the 
vortex core where the hydrodynamic theory does not hold. 



B. The superfiuid Magnus force 

In the superfiuid hydrodynamics one can refer this derivation to the superfiuid component with density p a . The 
Euler equation for the superfiuid component |^| after adding the external (5-function force F applied at the vortex 
line is 

Bv - F 

^ + (v s - VK = -V/i + -S 2 (f) , (10) 
dt p s 

where p is the chemical potential. 

For charged superfluids (superconductors) the Euler equation should include also the electromagnetic forces. In 
particular, the chemical potential must be replaced by the electrochemical potential. But except for vortex cores one 
may use the quasineutrality condition that the total electron charge is approximately equal to the background ion 
charge. This allows to neglect the chemical potential gradient. Finally the Euler equation may be written as 

Bv - e / -, 1 -> \ F 

+ (v s • VK = - (E+ -[*. x H] + -<5 2 (r) , (11) 
dt m \ c J Ps 

where E and H are the electric and the magnetic fields. 

Let us consider a vortex line in a neutral or charged superfluids with the velocity field Eq. ([!]). Now circulation is 
quantized, and the circulation quantum is k — h/m in a Bose superfiuid and n = h/2m in a Fermi superfiuid. One can 
repeat the analysis of the momentum balance for a cylindrical region around the vortex line. Now the velocity v is 
the superfiuid velocity v s far from the vortex line and the Bernoulli law is used for variation of the chemical potential 
or the electric potential in the neutral and the charged superfluids respectively. But the momentum of the superfiuid 
component is not conserved because of interactions with quasiparticles (and impurities in the case of superconductors) 
in the vicinity of the vortex line. One may assume that all these interactions are incorporated by the external force 
F, localized at the vortex line. Then one has instead of Eq. @: 

Ps[{v L -v s ) x Fz] =F . (12) 

The force F, which enters the theory as a (5-function force, is distributed over a small vicinity of the vortex line in 
reality. The dimension of this vicinity is not necessary to be of the vortex-core size, but must be smaller than all 
relevant hydrodynamic scales (e.g., the intervortex distance, or the vortex line curvature radius). 

Replacing the external force by the Magnus force, the Euler equations for the neutral and the charged superfluids 
are 

+ (v s ■ V)v s = -V/x + [{v L - v s ) x k]5 2 {t) , (13) 
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dv - e f - 1 - \ 

-g^ + {vs- V)v s — — yE + — [v s x B] J + [(vl — v s ) x K]S 2 (r f ) . (14) 

One can transform the Euler equation using the vector identity 

V)iJ s =V— -tf s x [X7xv s ] . (15) 

In a neutral superfluid vorticity is concentrated on the vortex line: [V x v s ] — K8{r). But in a superconductor 
V xv s = KS(r) - -^H. 

s v ' rnc 

Then the Euler equations are 

^■ = -^(f+f)+PiXi5]J 2 (f), (16) 



f = ^-v(f) (17) 

This analysis demonstrates that the total external force on the superfluid in the vicinity of the vortex line is exactly 
balanced by the superfluid Magnus force p s [{vL — v s ) x k]. In fact, the term oc vl in the Euler equation may be 
received from a pure kinematics: it presents the flow of the vortex lines across the line between two points which 
changes the phase difference between them (the phase slip). After replacing the external force by the Magnus force 
in the Euler equation, the latter does not contain any information on the nature and the magnitude of the external 
force. But the Euler equation is not sufficient for description of superfluid motion: an additional equation for the 
vortex velocity vl is necessary. In order to derive it, one should specify the force F. This may be done by considering 
the momentum balance of the whole liquid, but not only its superfluid component. 



C. Equation of vortex motion and effective Magnus force 



The momentum balance for the whole liquid in a cylindrical region around the vortex line has been studied in 
Ref. H using the collisionless kinetic equation for quasiparticles and assuming translational invariance. The balance 
yielded the equation of vortex motion which is a linear relation imposed on three velocities v s , v n , and vl- Since 
translational invariance was assumed, this equation depended only on velocity differences, but not their absolute 
values. However, one can also include into this balance some interactions with external fields, e.g., with impurities in 
superconductors. Assuming the axial symmetry in the plane normal to the vortex line, the momentum balance yields 
the following equation of vortex motion: 

Ps[(vl - v s ) X k] = -D(v l - v n ) -D'[zx (v L - v n )] - dv L - d'[z x v L ] . (18) 

Comparing it with Eq. (|l^), one obtains the expression for F: 

F = -D(v L - v n ) - D'[z x (v L - v n )} - dv L - d'[z x v L ] . (19) 

The forces proportional to D and D 1 are due to scattering of free quasiparticles by the vortex, therefore they are 
proportional to the difference between the drift velocity of quasiparticles (the normal velocity v n ) and the vortex 
velocity vl- The forces proportional to d and d! are due to interaction between the vortex line and impurities which 
are frozen into the crystal and therefore do not move if the crystal is at rest (The case when the crystal is not at rest 
is discussed in Ref. [ p6[ ) . These forces include also the interaction of impurities with the quasiparticles bound in the 
vortex core, and therefore moving with vl, but not with v n . 

One can rewrite the equation ( |l8| ) of vortex motion collecting together the terms proportional to the velocity vl'- 

Pm[vl x k] + rjv L = p s [v s x k}+ Dv„ + D'[z x v n ] . (20) 

Here pm = p s — (D' + d')/n and r/ = D + d. The forces on the right-hand side are driving forces produced by the 
superfluid and normal flows. In the theory of superconductivity the force Fl — p s [v s xk] = (1/c)[j s x $o], proportional 
to the superfluid velocity v s (or to the supercurrent j s = en s v s ), is called the Lorentz force. Here $o = hc/2e is the 
magnetic-flux quantum and the vector $o is parallel to k. There are also forces on the vortex produced by the normal 
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current j n = en n v n . One can find discussion of the effect of the normal-current force on electrodynamics of a type II 
superconductors in Ref. 

The left-hand side of Eq. (|2(]) presents the response of the vortex to these driving forces. The factor pm, which 
determines the amplitude of the effective Magnus force on the vortex is not equal to the superfluid density p s in 
general: it may be more or less than p s . Note that the Hall conductivity is governed by pm, but not by p s . At low 
magnetic fields the normal current is small compared to the supercurrent, i.e. the total current j ~ j s = en s v s and 



one may neglects the terms ex v n on the right-hand side of Eq. (20). On the other hand, the electric field is connected 
with the vortex velocity by the Josephson relation E = ~[H x Vl\- Then the equation of vortex motion is equivalent 
to the Ohm law connecting the current and the electric field. One can easily check that the Hall component of the 
conductivity is linear in pm- 

In the superfluidity theory they usually present the equation of vortex motion using the mutual friction parameters 
B and B' introduced by Hall and Vinen Because of translational invariance d = d! = for superfluids, and 
neglecting the normal motion (v n — 0) the equation is 

{■, Pn D A . Pn Dr „ PsPMK 2 -, . PsKD _ 

v L = [1- —B v s + —B[z x v s ] = -v s + —j— -[z x v s ] . 21 

V 2p J 2p D z + {p M K>y D 2 + {pMK) 2 

In the next section we shall calculate the amplitude D 1 of the Iordanskii force analyzing the interaction of the vortex 
with phonons in the Born approximation.. 



III. IORDANSKII FORCE AND AHARONOV-BOHM EFFECT 



A. Nonlinear Schrodinger equation and superfluid hydrodynamics 

The Gross-Pitaevskii theory |^9| has suggested the nonlinear Schrodinger equation to describe a weakly nonideal 
Bose gas: 

^=-^ V V + W^- (22) 

Here tp — aexp(i</>) is the condensate wave function and V is the amplitude of two-particle interaction. Using the 
Madelung transformation |47| , this equation for a complex function may be transformed into two real equations for 
the liquid density p — ma 2 and the liquid velocity v = ^rV^ where k — h/m is the circulation quantum. Far from 
the vortex line these equations are hydrodynamic equations for an ideal inviscid liquid: 

^+V(H = 0, (23) 



ov — -» 

— + {v-V)v = -Vp . (24) 

Here p = Va 2 /m is the chemical potential and c is the sound velocity. 

Suppose that a plane sound wave propagates in the liquid generating the phase variation <j>(r, t) = 4>o exp(ifc ■ r — iwt) . 
Then the liquid density and velocity are functions of the time t and the position vector r in the plane xy and can be 
written in the form 

p(f, t) = p + p (1) (r, t) , v(r, t)=v + %) (r, t) , (25) 

where po and vq are the constant density and velocity in the liquid without the sound wave, p(\)(f, t) and vt\)(f, t) — 

#-V<^> are the changes in the density and the velocity due to the sound wave. They should be determined from the 
hydrodynamic equations ( [23] ) and (|24|) after their linearization. In particular, Eq. ( p4[ ) gives the relation between the 
density variation and the phase 4>: 

Substitution of this expression into Eq. ( p^ ) yields the wave equation for a moving liquid. The sound wave has the 
spectrum lo = ck + k ■ vq. The sound propagation is accompanied with the transport of mass. In the reference frame 
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moving with the average velocity vo of the liquid the mass transport is determined by the mass current j ph which is 
of the second order with respect to the wave amplitude. Averaging over the wave period, one obtains: 



f h {k) = (p {1) v (1) ) = P^l^-k = n(p)p . (27) 



This expression supposes that the plane sound wave corresponds to a number n(p) of phonons with the momentum 
p = hk and the energy E = e(p) + p ■ vo where e(p) = cp is the energy in the reference frame moving with the liquid 
velocity vq. Then the total mass flow is 

j = Pom + tq / d 3 pn(p)p . (28) 



h 3 . 

In the thermal equilibrium at T > 0, the phonon numbers are given by the Planck distribution n(p) = no(E,v n ) 
with the drift velocity v n of quasiparticles: 

n (E 7 V n )= E(p)-p-v n 7 = e(p)+p-(v -v n ) 7* ( 29 ) 

exp j, — 1 exp w F T — ^ - 1 

Linearizing Eq. ( p9| ) with respect to the relative velocity vo — v n one can see that this expression is equivalent to 
the two-fluid expression j — pv s + p n (v n — v s ) = p s v s + p n v n assuming that p = po = p s + p n , vq = v s , and the normal 
density is given by the usual two-fluid expression: 

1 f dn a (e,0) 2 



3/i 3 J dE 

This simple analysis demonstrates that two-fluid hydrodynamics with phonon quasiparticles is identical to the 
nonlinear Schrodinger equation with thermally excited sound waves. A next step is to analyze scattering of phonons 
by the vortex in hydrodynamics of an ideal liquid. 



B. Scattering of phonons by the vortex in hydrodynamics 

The phonon scattering by a vortex line was studied beginning from the works by Pitaevskii ^] and Fetter [^9| . 
Interaction is weak and one may use the perturbation theory (the Born approximation). 

Let us consider a sound wave <j>(f, t) = tfio exp(ifc • r — iwt) propagating in the plane xy normal to a vortex line along 
the axis z. Then in linearized hydrodynamic equations of the previous section the fluid velocity vq should be replaced 
by the velocity v v (r, t) around the vortex line which now depends not only on the position vector r, but also on time, 
since the vortex line is not at rest when the sound wave propagates past the vortex. This means that r in Eq. (|l|) 
must be replaced by r — viX and dv v /dt — — ■ V)w„ = — V(ul • v v ). Since there is no external force on the liquid, 
the vortex moves with the velocity in the sound wave: vl = U(i)(0, t). In the presence of the vortex line linearized 
hydrodynamic equations are: 

^ + PoV • v (1) + v v ■ V P(1) = , (31) 



Since Va* = %^ P{x), Eq. (|3j) yields 



c%) ^ _ . dv v 

— = -V M -VK -v (1) ) - — 



(32) 



(33) 



Finally the linear equation for the phonon-induced phase is 



d 2 cj) 



C^V 



-2v v {r) ■ V 



dt 



<i>(r) - -m 



(34) 
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In the Born approximation one treats interaction with the vortex velocity field [the right-hand side of Eq. (|3J)] as 
a perturbation. Then 

4> = <j>o exp(— iujt) -|exp(ifc • r) + — J d,2r\H^ (k\f— fi\)k ■ v v (f\)[2 exp(ifc • r%) — 1] | (35) 

Here H^'(z) is the zero-order Hankel function of the first kind and jH^\k\f — r*i |) is the Green function for the 2D 
wave equation, i. e. satisfies to the equation 

(fc 2 - V 2 )0(r) = tf 2 (f - n) ■ (36) 

In our problem the small perturbation parameter is nk/c which is on the order of the ratio of the wavelength 2ir/k to 
the vortex core radius r c ~ n/c. 

The standard procedure in the scattering theory is the following. One uses the asymptotic expression for the Hankel 
function at large values of the argument: 

lim HP(z) = «/X e i(*-*/4) . (37) 

Then it is assumed that the perturbation is confined to a finite vicinity of the line, where r\ <C r, and 

\f-fi\ «r- . (38) 

r 

Finally the phase field at large values of kr may be presented as a superposition of the plane wave oc exp(zfc • r) and 
the scattered wave oc exp(ifcr): 



4> — 4>q exp(-iujt) 



expufc ■ r) H -=- exp(ifcr) 

\/r 



(39) 



Here a((p) is the scattering amplitude which is a function of the angle ip between the initial wave vector k and the 
wave vector k! — kf/r after scattering. For scattering of phonons by the vortex the Born approximation yields that 
@ 

, . [k 1 j7r rA ? n j* ( 1 1 \ 1 [k k ^smipcostp , An . 

a(ip) = \ — e 4 \k x k'\ ■ k - —r = -\ e' 4 , (40) 

^' \ 2n c J \q 2 2k 2 J 2 V 2tt c 1 - cos tp ' v ' 

where q = k — k' is the momentum transferred by the scattered phonon to the vortex, and q = 2fcsin(<y9/2). The 
second therm I /2k 2 in parenthesis is due to the vortex line motion |50| . 

In order to find a force on the vortex from the sound wave, one must determine a contribution of the wave to the 
average total momentum flux Ff h — J dSjl^j over the cylindrical surface around the vortex line where 

nf/ = (P(2))5ij + (p(l)(v {1) )i)v VJ + (P(l)(W(l))j>m + Po{(u(l))i(W(l))j) • (41) 

2 

Here P( 2 ) = Pof i (2) + §^~2^ ^ s *-^ e second-order phonon contribution to the pressure where = po/c 2 . According to 

2 

the Euler equation (Q) the second-order contribution to the chemical potential is /X(2) — — Then 

In Appendix ^ it is shown that if the perturbation of the sound wave is confined to a finite vicinity of the line, the 
force on the line from the sound wave, 

FP h = <J llC jP h -a ± c[zxf h } , (43) 
is determined by two effective cross-sections: the transport cross-section for the dissipative force component, 



°ll = / (7 (V 3 )( 1 - cos<p)dip , (44) 
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and the transverse cross-section for the transverse force component, 

a± = J a(f) sin ip dip . (45) 

The differential cross-section cr((p) = | a(<^>) | 2 in these expressions is known due to Eq. ( [fb"| ) for the scattering amplitude 
a(ip). It is quite natural that in the Born approximation the transverse cross-section vanishes since the differential 
cross-section is quadratic in the circulation re. 

However, the standard scattering-theory approach fails to describe the phonon scattering at small angles ip. Indeed, 
the velocity v v induced around the vortex is decreasing very slowly, as 1/r. Therefore the terms oc v v in the phonon 
momentum flux, Eq. (|4l]), are important in the momentum balance. In addition, the scattering amplitude is divergent 
at ip — > 0: 



lim a((f) = 




This divergence is integrable in the integral for the transport cross-section, Eq. (Hj). So the calculation of the 
transport cross-section is reliable. Contrary to it, the integrand in Eq. (^5|) for the transverse cross-section has a pole 
at if = 0, and the contribution of this pole requires an additional analysis. A proper analysis of the phonon small-angle 
scattering was fulfilled in Refs. |||^J. However, in the recent publication Demircan, Ao, and Niu |52| considered the 
phonon scattering by the vortex ignoring special features of the small-angle scattering. This is the reason why they 
could not find the transverse force from phonons on the vortex. 



C. Small-angle phonon scattering in the Born approximation and the Iordanskii force 



At small scattering angles <p < 1/y/kr the asymptotic expansion given by Eq. ( |39| ) does not hold. The accurate 

for small angles was done in Ref . || . A simplified version of this calculation is 
;hat at ip -C 1 



calculation of the integral in Eq. 



presented in Appendix [Bj. It yields t 

(f) = c/)q cxp(—iujt + ik ■ r) 
Using an asymptotic expression for the error integral 

z 

2 



d>(z) 



e _t dt 




oo 



(47) 



exp(— z 2 



(48) 



at large z, one obtains for angles 1 ^> ip ^> l/Vkr 



4> = 0o ex p(~ iwt) 



i-T ( f \ ik k 1 ( 

eMlk • r K 1 + ^Rj + ~V^ exp r 



IT 

tr + i- 
4 



(49) 



The second term in brackets coincides with scattering wave at small angles ip <C 1 when the scattering amplitude is 
given byEq. @). But now one can see that the standard scattering t heory misses to reveal a very important non- 
analytical correction to the incidental plane wave. We shall see in Sec. Ill E that the factor ±nk/c which determines 
this correction is exactly the phase shift of the sound wave along the quasiclassical trajectories past the vortex on the 
right and left sides. This is a manifestation of the Aharonov-Bohm effect JjJ: the sound wave after its interaction with 
the vortex velocity field has different phases on the left and on the right of the vortex line, and this phase difference 
results in an interference. 

In the interference region the velocity induced by the sound wave is obtained by taking the gradient of the phase 
given by Eq. ([47|). The velocity component normal to the wave vector k, 



n dip . k k 

; (l)-L = n — "a = W exp(-iu;t + ikr)- — 



(50) 



determines the interference contribution to the transverse force: 
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/ 



dSjpo{(v( l ))±{v( 1 )) j ) = j p ((v^)x(v(i))r)rd(p = ^^ Po4>l— —Vkr j dip cos l-krtp 2 ) = -^Po$—^- ■ (51) 



Note, that this force contribution arises from the interference region with the transverse dimension d int ~ y/rg/k. Here 
ro is the large distance from the vortex line where the momentum balance is considered. But the interference region 
corresponds to very small scattering angles ~ di nt /ro = 1/ ykr~o. Thus an infinitesimally small angle interval yields 
a finite contribution to the transverse force. One couldn't reveal such a contribution from the standard scattering 
theory using the differential cross-section. 

Exactly the same contribution to the transverse force, as in Eq. (|5l|), arises from the term v Vi {p^(v^)j) in 
the momentum-flux tensor, Eq. (^|). In this term the mass flow (p(i)(u(i))i) for the plane wave may be used [see 
Eq. (p7|)], since we take into account only terms which are of the first-order in the Born parameter nk/c. Finally this 
yields the transverse force given by the transverse cross-section cti = k/c which is linear in the circulation quantum 
k and therefore cannot be obtained from the differential cross-section quadratic in n. This is the Iordanskii force. We 
shall rederive it in the next section using expansion in partial waves in order to demonstrate the analogy with the 
Aharonov-Bohm effect. 



D. Partial- waves analysis and the Aharonov-Bohm effect 

Neglecting the terms of the second order in v v , the equation 

k 2 <j) ~ f-iV + -v v \ = (52) 



c 



is equivalent to Eq. ( p4[) written for the harmonic sound wave with frequency u) = ck. The contribution from 
the vortex-line motion [the term cx 0(0) on the right-hand side of Eq. (p4|)] is now neglected as unimportant for 
the transverse force. Equation (p2|) is analogous to an equation which describes interaction of an electron with the 
magnetic flux <E> confined to a thin tube (the Aharonov-Bohm effect |Q): 



Eip(r) = — -ihV - - A Mr) . (53) 

Here ip is the electron wave function with energy E and the electromagnetic vector is connected with the magnetic 
flux by the relation similar to that for the velocity v v around the vortex line [Eq. (|l])]: 

(54, 

Let us look for a solution of this equation as a superposition of the partial cylindrical waves using the cylindrical 
system of coordinates (r, ip) : 

ip = y]i/>i(r) exp(ilip) . (55) 
/ 

The partial-wave amplitudes should satisfy equations 

~r~o~ H 1 5 — + k ipi = . (56) 

dr z r ar r z 

Here k is the wave number of the electron far from the vortex so that E = h 2 k 2 /2m and 7 = $/3>i where $1 = hc/e 
is the magnetic-flux quantum for one electron (two times larger than the magnetic-flux quantum $0 = hc/2e for a 
Cooper pair). A solution of this equation is the Bessel function J|;_ 7 |(fcr) with the following asymptotics at large 
arguments: 



J |; _ 7| (fcr)^y -J-cos(fcr-||Z- 7 |-|) . (57) 
On the other hand, the expansion of the plane wave in the partial cylindrical waves is 

exp(ik ■ r) = exp(ikr cos if) = Ji(kr) exp[il(ip + 7f/2)] , (58) 
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or at large kr: 



expfzfc ■ r) = \ cos (kr I \ exp[il((p + ir/2)] . (59) 

V irkr V 2 4/ 



In order to obtain the solution "the incoming plane wave + the scattered cylindrical wave" like Eq. (|39|), one should 
determine the partial waves from the condition that the incoming components cx exp(— ikr) in the plane wave and 
in the solution Eq. ([55|) coincide. This yields that 



.7T 



2 



7T , , 7T 



- \ l - l\) cos(fcr - -\l - i\ - -) . (60) 



2' " 4 



Then the solution of the Aharanov-Bohm problem takes the form of Eq. (p9h with the scattering amplitude 



= V 2^fc 6XP ~ exp ( 2i(Si )] ex P(*M , (61) 

7T ; 

where 

S t = (/ - ]/ - 7|)tt/2 (62) 

is the partial- wave phase shift. 

Equation ( |45| ) for the transverse cross-section may be rewritten as an expansion in partial waves f3l| : 

c_l = / |a((p)| 2 sin^rfy 3 — r J]] sin2(<5; — <5;+i) ■ (63) 



Using the phase shift values for the Aharonov-Bohm effect, Eq. (|62j), the transverse cross-section is 

cr_L = sin 27T7 . (64) 

k 

One can obtain the cross-section for phonon scattering from this expression assuming that 7 = —nk/2irc and expanding 
the sine- function in small 7. 

The cross-section for the Aharonov-Bohm effect is periodical in the magnetic flux with the period equal to the 
one-electron flux quantum <!>i. If the electron is scattered by the Cooper-pair flux quantum $0 = $1/2 the transverse 
cross-section vanishes. But presented analysis of the Aharonov-Bohm effect is based on the assumption that the total 
magnetic flux is concentrated in a very thin tube. Namely, the radius of the tube must be much less than the electron 
wavelength. This condition certainly doesn't hold in superconductors, where the electron wavelength is on the order of 
the interatomic distance and the effective radius of the magnetic-flux tube is the London penetration depth. Therefore 
one should consider scattering of electrons by a thick magnetic- flux tube f5l)| . Scattering of the BCS quasiparticles 
by the magnetic field of the vortex is contributes to the transverse force on the vortex, but this contribution cancels 
with the bulk electromagnetic force on the superfluid electrons || . 

In their original paper |jj Aharonov and Bohm considered only the effect of the magnetic-flux tube on the electron 
wave. The force from the electrons on the fluxon was considered later by Aharonov and Casher f52]| and therefore 
is called the Aharonov-Casher effect. So for our problem of interaction between quasiparticles and the vortex both 
effects, Aharonov-Bohm and Aharonov-Casher, are important. But these effects are so close one to another that 
throughout the paper we use only the name "Aharonov-Bohm effect" . 



E. The transverse force in the quasiclassical theory 



Now we shall show that the Iordanskii force follows also from the quasiclassical theory of scattering by the vortex, 
despite one cannot use the quasiclassical theory for phonons. But the quasiclassical theory is valid to describe the 
roton contribution to the transverse force. 

Let us consider a quasiparticle with an arbitrary spectrum s(p). If the quasiparticle moves in the velocity field 
induced by the vortex, its energy is E(p) — s(p) +p ■ v v which may be treated as a Hamiltonian to write the classic 
equation of motion: 
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dr BE p 

-77 = -q^ = vg- +v v , 

at op p 

% = -w = -^'*° )i (65) 

where vg{p) = de/dp is the quasiparticle group velocity. As well as in phonon scattering, we assume that the 
roton moves in the plane xy normal to the vortex. From this equations one can find the classical trajectory of the 
quasiparticle moving past the vortex line. Usually it is close to a straight line, and a distance of the straight trajectory 
from the vortex line is the impact parameter b. We are looking for the transverse force only, then we need only the 
variation Sp± of the momentum component normal to the initial momentum p which results from quasiparticle motion 



past the vortex. In the first order with respect to v v Eqs. (65) yield that 



— oo 

where I is the coordinate along the trajectory. On the other hand, the momentum of the quasiparticle is connected 
to the classical action: p — dS/dr. Then Sp±(b) = dSS(b)/db where the total variation of the classical action along 
the trajectory is a function of the impact parameter b: 

oo 

— OO 

The scattering angle of the quasiparticle is (p ~ —Sp±(b)/p and the transverse cross-section is 

oo oo oo oo 

f „ • / „ „ v SS(-oo) - SS(+oo) 1 f „ d f „ Kb k 



For phonons the group velocity vg is equal to the sound velocity and Eq. ( p8[ ) yields the correct value of a± even though 
the quasiclassical theory is not valid for phonons: there is no well-defined classical trajectories for phonons except for 
large impact parameters b at which the scattering angle ip is negligible. But in order to obtain the transverse cross- 
section <7j_, one should know only the action variation at large impact parameters whereas the derivative d5S(b)/db 
of the action variation, which determines the scattering angle, is not essential. 

It is important to note that the double integral of Eq. ( |6"§| ) is improper: its value depends on what integration 
is done first. The correct procedure which was justified in Ref. fj is to integrate along the trajectory first, and to 
integrate over the impact parameters afterwards. A way to check it is the following. We choose some finite limits in 
the double integral of Eq. ( |38| ) which means that the integration is restricted by some area around the vortex line. 
The integral depends on the shape of this area. For example, a circular border of the area yields a± by a factor 2 less 
than that of Eq. (^8|). But the full solution of the collisionless kinetic equation for the quasiparticles made in Ref. || 
showed that other terms contribute to the momentum balance, also originating from the slow decrease of the velocity 
field. Taking into account all of them, we arrive again to the expression for the force via the cross-section given by 
Eq. (|8]). The order of integrations in Eq. (^38j) assumes that the integration area has a shape of a rectangular with 
a long side along the quasiparticle trajectory. For such a shape all other contributions to the transverse force exactly 
cancel. 

The transverse cross-section crj_ determines the amplitude D' of the transverse force on the vortex in Eqs. ( |l9| ) and 
( pp| ) . For quasiparticles moving at an arbitrary angle to the vortex line v G in Eq, ( |68| ) is the component of the group 
velocity in the plane normal to the vortex line. The amplitude D' of the transverse force must be determined using the 
Planck distribution for phonons with the drift velocity v n . This distribution should be linearized with respect to the 
relative velocity v s — v n = (v s — v£) + (vl — v n ). But the velocity v s — vl enters the energy E{p) = e{p) +p- (v s — vl) 
of the quasiparticle in the reference frame connected with the vortex. Using the condition of the detailed balance one 
can show that the quasiparticle distribution function which depends only on the energy cannot produce a force. Thus 
the only part of the distribution function which contributes to the force is that linear in the relative drift velocity 
Vn ~ vl ■ Finally the amplitude of the transverse force is 

, 1 f dn o (e,0) 2 _ 
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Using the cross-section a± from Eq. (pq) . one obtains that D' = —np n . This means that in the translationally 
invariant superfluid without scattering by impurities the amplitude pm of the effective Magnus force on the vortex 
in Eq. (|2^) is equal to the total density p, but not to the superfluid density p s . Then neglecting the viscous forces 
[D = d = d! = in Eq. (|2C|)1 the vortex moves with the center-of-mass velocity. 

A rather simple and universal expression D' — —np n for the Iordanskii force amplitude tempts to claim its universal 
topological origin, since k in this expression is a topological charge. However, in the next section we shall see that 
the expression is not universal, in fact. For quasiparticles in a BCS superconductor with energy much exceeding the 
gap an additional small factor should be put in this expression. 

IV. IORDANSKII FORCE FOR QUASIPARTICLES IN BCS SUPERCONDUCTORS 

The wave-function of quasiparticles in the BCS theory has two components, 

^=(;$) . (70) 

which are determined from the Bogolyubov-de Gennes equations: 

fi 2 

- 77- f V 2 + kl) u(r) + A expUemvtf) = Eu(r) , (71) 
2m V / 



2m 



(v 2 + k%) v(r) + Aexp[-i6(r)}u(r) = Ev(r) . (72) 



Here Uf is the Fermi wave number. We neglect the magnetic field effect which is weak if the London penetration 
depth is large compared to other relevant scales. Without the vortex the order parameter phase 6 is a constant and 
the equations yield the well-known BCS quasiparticle spectrum E = \J '£ 2 + A 2 , where £ = (ft /2m) (k 2 — k 2 F ) is the 
quasiparticle energy in the normal Fermi-liquid. 

In Refs. |^,[n^] the Bogolyubov-de Gennes equations for a quasiparticle passing a vortex were solved with help of 
the partial-wave expansion, earlier used also in Ref. [^l). Quasiparticles with the energy close to the gap (£ <C A) 
behave as rotons and the transverse cross-section for them is given by Eq. ( p8| ) in which the group velocity for the 
BCS quasiparticles is vq = vf^/E. Here vf — hkp/m is the Fermi velocity. In the case when the quasiparticle 
energy is much more than the superconducting gap, the theory yielded that the transverse cross-section differed from 
the quasiclassical result of Eq. (pS) by the factor A 2 /2£ 2 . Now we rederive this result using the Born approximation 



similar to that for phonons in Sees. Ill B and III C 



We use the perturbation theory with respect to the gap A and the gradient of the order parameter phase V#. 
Then in the zero-order approximation u = uoexp(ifc • r) and v — 0. In the first-order approximation the second 



Bogolyubov-de Gennes equation (72) yields 



f Aexp(-i6>) Aexp(-i60 ft 2 .- - s l , r 

The first term in curled brackets yields a correction to the quasiparticle energy cx A 2 , but does not contributes to 
scattering which is determined by the order-parameter phase gradients. So we keep only the second term proportional 
to W. Inserting it to the first Bogolyubov-de Gennes equation ((n]) one obtains the following equation for the 
first-order correction to the plane wave: 



A 2 

(V 2 - k 2 )u {1) - (k ■ V6)—uv exp(ik ■ r) . (74) 



This equation is similar to equation ( |34| ) for the sound wave and using this analogy one easily obtains the expression 
for the transverse cross-section: 

A 2 7T A 2 K 

2£ z Kf 2£ z vf 

Since the group velocity of quasiparticles with E 3> A is about the Fermi velocity vf, this expression differs from 
Eq. ( pq ) by the factor A 2 /2£ 2 . Thus the Iordanskii force becomes small close to T c : it decreases proportionally to 
p s ~ A 2 , like the superfluid Magnus force. 
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V. MAGNUS FORCE IN THE JOSEPHSON JUNCTION ARRAY 



In the continuum limit, the equation of motion for a vortex has been derived in Ref. fl36| . This derivation has not 
revealed any force normal to the vortex velocity. Absence of the Magnus force suggests that the vortices move parallel 
to the driving force, i.e., normally to the current, and there is no Hall effect. Then in the limit of weak dissipation 
the ballistic vortex motion is possible, which is a free vortex motion without friction and the driving force. The 
Hall effect and the ballistic motion are incompatible, since the latter is a regime with a finite electrical field and no 
external current which is impossible for a finite Hall resistance. Though there have been experimental evidences of 
the ballistic vortex motion p4|| , one may suspect that a more sophisticated theory would reveal the Hall resistance, 
however small. In the present section it will be shown that the Hall effect is exactly absent in the classical limit for 
the J J A. It directly follows from the symmetry of the dynamic equations. 

Let us consider a conductor in a magnetic field H. When its symmetry is not less than the three-fold (which 
includes a triangular and square lattices), the Ohm law is: 

E = p L I + p H nx I , (76) 

where ft = H/H , E is the electrical field, pz, is the longitudinal resistance and pn is the Hall resistance in the magnetic 
field H . Now let us consider the transformation in which the directions of the fields E and H and the current / are 
reversed: 

E->-E, ft — ► —ft(H — > —H), f ->■-!. (77) 



The Ohm law Eq. (76) is invariant with respect to this field- current inversion only for a system without the Hall 
effect (ph = 0). On the microscopical level the field-current-inversion invariance is a direct result of the particle-hole 
symmetry which was shown to forbid the Hall effect in the Ginzburg-Landau theory (see |2^] and the references 
therein). 

Next we consider the JJA with the energy 

£ = ^E[^S^-^ sin %-^)] - (78) 



l.k 



and the equations of motion 



£ c w~sr - lc £ 8in <*<- - *<+*> + £ "o+A+i. = ■ < so > 

p p p 

Here VJ-is the electric potential, Qf = J2p ^TT+p^T+p ^ s ^ nc electric charge, and (f>f{t) is the gauge invariant phase at the 
node specified by the discrete two-dimensional position vector I, Ej is the Josephson coupling energy, lc = 2eEj/h is 
the critical current, C ; -- and ctj-- are the capacity and the conductance matrices respectively. In the external magnetic 

field H — V x A, the gauge invariant phase <fif is not single-valued; in fact only its difference between neighboring 
nodes is well-defined: 

2tt /•('+#) _ _ 

<t>T + p - <t>T = n + p 4 A - dl . (81) 

Here ip^ is the canonical phase at the node I, and the integral over A is taken between the centers of the two neighboring 
nodes. The canonical phase is not single- valued too, but its circulation along any closed path through the nodes of JJA 
is always an integer number of 2-7T, while the circulation of the gauge invariant phase may be any number depending 
on the magnetic field. 

When both the external current I and the magnetic field H are applied to the JJA, the gauge-invariant phase 

can be presented as d>f = d>S + d>l. Here iff is the time- independent phase in the equilibrium state without an 
L i i i 

external current, and <f>l is the time-dependent contribution to the phase from the external current /. The multi- 
valuedness of the phase related to the magnetic field is present only in the static phase <f)S: the time-dependent 
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dynamical contribution <pl to the phase is single-valued. One sees then that the field-current inversion [Eq. ( |77 
simply corresponds to the change of signs of all phases, and the equations of motion, Eqs. ( f79f ) and (|so|), are invariant 
with respect to this transformation. It proves that the Hall effect does not exist in the JJA, i.e., the effective Magnus 
force vanishes. 

The crucial point of this very simple derivation is that one can use the static vortex solution for the dynamical 
problem. This assumes that singularities of the phase distribution related to the presence of vortices are kept at 
rest despite the vortices themselves are driven by the Lorentz force. For continuous superconductors this approach 
is invalid and our derivation does not work. So there is a fundamental difference between vortices in a lattice and 
vortices in a continuous superconductor. Indeed, in the lattice there are no singular vortex lines. They appear only 
in the continuum limit. At best, one can define the lattice cell containing the vortex center. This definition has 
been borrowed from the continuous theory: it is the cell, around which the circulation of the phase ip is equal to 
2tt. However, in the lattice the circulation around a closed path is not well-defined. Let us consider some closed 
path through a discrete number of nodes with the phase circulation 2tt. One may change the phase difference by — 2n 
between any two neighboring nodes on the path without any effect on observed physical parameters (currents, voltages 
and so on). Then the circulation vanishes along the path considered, but must appear along a path over other nodes. 
Thus one cannot locate the position of the phase singularity. In order to avoid this ambiguity in the JJA model, a 
special rule has been formulated: the phase difference between two neighboring nodes must not exceed it. When for 
some bond the phase difference achieves the value tt, one must redefine the phases; as a result, the vortex center is 
put into another cell. This procedure is usual for numerical studies of the vortex motion in JJA |39f| . However, this 
rule is not obligatory for the dynamic theory of JJA. Instead, one may keep 2ir circulations of the phase ip at fixed 
cells during the dynamic process without worrying where the vortex center (defined according to the aforementioned 
rule) is really located. 

It is important to stress, that one cannot derive the effective Magnus force in JJA using the continuous approach. 
Let us discuss this in more details. In the continuum limit the set of discrete vectors I is replaced with the continuum 
space of I. One can define the field of the canonical (but not gauge invariant!) phase (f(l) in this space everywhere 
except for the singular points which are the centers of the vortices with the phase circulation 2tt around them. It 
is assumed that the spatial variation of the phases is small and the phase differences can be replaced by the phase 
gradients according to 

( 82 ) 

Then one can derive the partial differential equations for a continuous field of the canonical phase f(l,t) which 
correspond to some Lagrangian L{(p(l,t)} (the dissipation is neglected now). The Lagrangian may include the term 
proportional to the time derivative of <p(l,t) which is called the Wess-Zumino term ]20| : 

L{<p(T, *)} = + LoW(l, t)} . (83) 



It is possible to derive the equation of vortex motion from this Lagrangian following Refs. J3 6|j3 7] ] . One must use the 
phase field for a slowly moving vortex: 

^(M) = arctan^^, (84) 
l x - x(t) 

where r(t) — [x(t),y(t)] is the 2D position vector of the vortex center. Substituting the vortex solution into the 
Lagrangian density given by Eq. ( p3|) and integrating over the xy plane, one obtains the effective Lagrangian which 
is now a functional of the trajectory for a moving vortex: 

C v {r{t)} = -Trgf • [z x r\ + C% {r(t)} . (85) 

Varying this Lagrangian with respect to r(t), one obtains the equation of vortex motion with the effective Magnus 
force oc q: 

2nq[?x z] = F s . (86) 

Here Fs includes all other forces (the Lorentz force and the inertia force) obtained from the Lagrangian Cq {f(t)} 
without the Wess-Zumino term. However, the factor q is not defined. If q is constant, it has no effect on the field 
equation for ip(l, t) since the Wess-Zumino term is a full time derivative in the field Lagrangian, Eq. (jS3|). The unknown 
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factor q should be proportional to some electric charge, since the charge is a variable conjugate to the canonical phase 
ip. But it remains unclear what is this charge: either the background charge determined by the whole Fermi-see of 
the superconducting island, or an external charge induced outside as suggested in Ref. [Q. Thus in the continuum 
limit the problem of the effective Magnus force and the Hall effect in the JJA remains unresolved. It must not be a 
surprise since in the continuum limit the JJA model becomes translationally invariant and "forgets" that originally 
it had been a lattice model without translational invariance. Meanwhile, the latter is crucial for the amplitude of the 
effective Magnus force. An additional physical principle beyond the continuum theory should be involved to obtain 
the equation for the vortex velocity. This principle is provided by our symmetry analysis. According to it q — and 
one should not include the Wess-Zumino term into the field Lagrangian. 

The presence of the external charge has no effect on our symmetry analysis. In order to take into account the 
external charge, one should use the Gibbs potential G = £ — V ex Y^fQfi where £ is given by Eq. ( f78| ) and V ex is 
the electric potential which creates the external charge Q ex — V ex J2p ^T+a- Then introducing the charge deviation 
Q'r = Qf~ Q ex -i one returns back to the energy £ with Q ~ instead of Qf. These arguments show that the external 
charge cannot lead to the Hall effect: its effect is restricted with the shift of the Fermi level, but the particle-hole 
symmetry is restored with respect to the new Fermi level. However, the external electric charge may produce the 
Magnus force in the quantum theory of JJA which takes into account the electron charge quantization . Then the 
Magnus force and the Hall conductivity are periodic in the electron charge. 

VI. SUMMARY AND DISCUSSION 

We have shown how the Magnus force appears in the equation of motion for a superfluid component (the superfluid 
Magnus force) and the equation of motion for a vortex (the effective Magnus force). Whereas the superfluid Magnus 
force proportional to the superfluid density is known exactly (from classical hydrodynamics, or from the Berry phase 
approach), there is no general expression for the effective Magnus force: it depends on interaction of the vortex with 
quasiparticles and with the external fields, like those from impurities in a dirty superconductors. Meanwhile, it is 
mostly the effective Magnus force which determines the observable effects: the mutual friction in superfluids, the Hall 
effect and the acoustic Faraday effect in superconductors, vortex quantum tunnelling. 

We have calculated the contribution of quasiparticles to the effective Magnus force for phonons in a superfluid 
and for high-energy quasiparticles in a BCS superconductor using the Born approximation. The transverse force 
from quasiparticles on the vortex (the Iordanskii force) originates from interference between quasiparticles passing on 
different sides of the vortex (the Aharonov-Bohm effect). 

Our symmetry analysis of the classical Josephson junction array has demonstrated that the effective Magnus force 
exactly vanishes and there is no Hall effect despite the finite superfluid density. One may formally interpret this result 
that the force from external fields violating translational invariance exactly compensates the superfluid Magnus force, 
though the analysis is not able to reveal these two forces separately, but only their joint outcome, namely, the effective 
Magnus force. 

The Ao-Thouless theory yields only the superfluid Magnus force which appears in the momentum balance of the 
superfluid component (the condensate). Indeed, Gaitan |2C]] derived the Ao-Thouless result for a charged superfluid, 
analyzing the momentum balance for the condensate. In order to derive the effective Magnus force (the total transverse 
force on the vortex), one must consider the momentum balance for the whole system. 

Ao, Niu, and Thouless J53| stated that the Iordanskii force did not appear in their Berry phase approach. But the 
recent paper on the Born quasiparticle scattering by Demircan, Ao, and Niu |?2[ demonstrated that they admitted 
the force from quasiparticle scattering in their approach, but concluded that this didn't yield any transverse force on 
the vortex. This conclusion was based on a wrong analysis of the Born phonon scattering missing the contribution 
from the Aharonov-Bohm interference. Thus the source of controversy is not in a difference of approaches, but in the 
problem how to calculate integrals for the Born phonon scattering. 

The Ao-Thouless theory rejects also any force on the vortex from the external fields, like the Kopnin-Kravtsov force 
in a dirty superconductor. On the basis of this theory Gaitan and Shenoy |2^] predicted the finite effective Magnus 
force and the Hall effect for the Josephson-junction array. This prediction contradicts to our symmetry analysis and 
to the experiment. Gaitan and Shenoy |2|] used in their analysis the Wess-Zumino term in the Lagrangian for the 
continuum limit of JJA. We have shown in Sec. why this approach is not reliable. On the other hand, Zhu, Tan, and 
Ao |Q try to conciliate the negligible Hall effect in JJA with the Ao-Thouless theory supposing that the vortices cross 
JJA between superconducting islands where the superfluid order parameter vanishes and therefore the Magnus force 
is negligible. But the superfluid order parameter does not vanish between the islands even though it is much smaller 
than in the islands. In fact, this small order parameter determines a finite superfluid density in the continuum limit. 
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Thus in JJA the Hall effect vanishes despite a finite superfluid density, contrary to the prediction of the Ao-Thouless 
theory. 

Makhlin and Volovik p5| suggested that the superfluid Magnus force in JJA is nearly compensated by the force from 
the bound states in the junctions (the spectral flow of bound states). But they did not conclude that the compensation 
is complete, and assumed the Fermi superfluid in islands and the SNS Josephson junctions. Our analysis shows that 
the Magnus force exactly vanishes for JJA independently on microscopic nature of the superconducting islands and the 
junctions. This shows that the core bound states and the spectral flow are not the only explanation for compensation 
of the Magnus force in the systems without translational invariance. 
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APPENDIX A: THE FORCE ON THE LINE SCATTERING THE SOUND WAVE 



First we derive the analogue of the optical theorem for the sound wave. For the latter we use the asymptotic 
representation, Eq. (39), in which the scattering amplitude a(ip) is not necessarily obtained in the Born approximation. 
But in general a(p>) should satisfy the condition that the total mass flow through the cylindrical surface surrounding 
the scattering line vanishes. 

An asymptotic expression for the average mass flow from the sound wave is 



J 



(P(1)V 



(1)) = PoPo 



K 2 k 

8^c 



- \ar - - - 1 

k + —k' - (k + fc')^[Im{a} cos(fcr - k ■ f)] + Re{a} sin(fcr - k ■ f)] 
r Jr 



(A.l) 



The condition that the total flow through the cylindrical surface around the scattering line vanishes is 

K 2 fc 2 f f . \ a(<p)\ 2 Im{ a (yQ} 



jf'dSi = I (p^v Wi )dSi = p <f>l 



8tt 2 c 



cos <p 



- (1 + cos ip) cos[fcr(l — cos tp)] > dp = 



(A.2) 



The integral over the term oc Im{a(p)} expands only over the region of small angles since kr 3> 1. Finally this 
condition yields 



- 2^lm{a(0)} + J\a 



(p)\ 2 dp = Q 



(A.3) 



Next let us consider the momentum balance which determines the force on the scattering line from the sound wave 
F? h = -fdSjn% where: 



n 



ph 



? (ph) 



Po- 



Sij +Po((«(i))i(«(i))i) 



(A.4) 



\Po 2 2 t 

The pressure term vanishes after averaging, but the convection term is essential and yields for the force on the vortex: 



K 2 k 
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k C0S v + M^ fcr -_ Im H°» cos ( l w , 2fc 
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r ' V k 



(A.5) 



With help of the optical theorem Eq. ( |A.3| ) one obtains the expression Eq. (|4J) with the effective cross-sections 
determined by Eqs. (E|) and (E|). 
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APPENDIX B: SMALL-ANGLE SCATTERED SOUND WAVE 



Using the asymptotics of the Hankel function, Eq. (pq) can be rewritten as 



o exp(— iwt) I exp(ifc • r) H \ - — / d%?\ exp(i/c • r\ + ik\r — 7 r i I) 

c V 2nr 



.. k ■ [z X Fi] 



(B.l) 



Here the effect of the vortex- line motion was neglected as irrelevant for small-angle scattering. Expansion Eq. (38) is 
not accurate enough and next terms of the expansion must be kept: 



r — ri i=a r 



(fi ■ r) r\ (f\ ■ t)* 



2r 2r 3 



(B.2) 



The terms of the second order in r\ are important since the perturbation is not well localized near the vortex line, 
but decreasing slowly when n is increasing. Using the Cartesian coordinates of the position vector r\(x, y) and the 
inequality ip <C 1, one obtains 



i exp(— iwt) < exp(ifc • r) + 



nk 2 



2nr 



dx dy exp 



ik [ r — yip + 



y 



x 2 + y J 



(B.3) 



The double integral in this expression may be transformed into the error integral: 



oo oo 



dx / dy exp 



ik I r — y(f + 
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Then Eq. (B.3) coinsides with M 
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